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Abstract. There are several known schemes for entangling trapped ion quantum bits 
for large-scale quantum computation. Most are based on an interaction between the 
ions and external optical fields, coupling internal qubit states of trapped-ions to their 
Coulomb-coupled motion. In this paper, we examine the sensitivity of these motional 
gate schemes to phase fluctuations introduced through noisy external control fields, 
and suggest techniques to suppress the resulting phase decoherence. 
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Scalable quantum computing presents a direct application for the study and control 
of large-scale quantum systems. The generally accepted requirements for quantum 
hardware pQ include identifiable two-level-systems for storing information in the form 
of quantum bits (or qubits), and methods to externally manipulate and entangle qubits 
through quantum logic gate operations. The implicit interconnects represented by 
entangled quantum systems provide the power behind quantum computation, giving 
rise to certain applications such as Shor's factoring algorithm [2] and Grover's search 
algorithm [Hj that exceed the capabilities of classical computers. However, in engineering 
complex entangled states of many qubits, it is critical to control the phase of the system 
of qubits and the phase of the classical control parameters that guide the quantum gates. 
Techniques of quantum error-correction [3J |3| appear essential to stabilize quantum 
computations, but to reach fault-tolerant error-correction thresholds [Hj, the host system 
must already possess a great deal of passive stability and must be relatively insensitive 
to external noise. 

One of the most promising quantum computing architectures is a system of cold 
atomic ions confined in free space with electromagnetic fields jTJ |HJ El- Here, qubits 
are stored in stable electronic states of each ion, and quantum information can be 
processed and transferred through the collective quantized motion of the Coulomb- 
coupled ion crystal. Applied electromagnetic fields (usually from a laser) enable this 
coupling between internal qubit states and external motional states, following several 
known quantum gate schemes [TUl HH E2]- Trapped ion quantum gates are thus 
highly susceptible to noise on the applied laser fields in addition to ambient electric and 
magnetic fields. In most cases, the relevant phases of the laser fields for quantum gates 
can be sufficiently stable during the evolution of a given gate, with gate speeds typically 
faster than about lOO/xs. However, for extended operations involving many successive 
gates, it will be difficult to maintain optical phase stability over the duration of the 
quantum computation. 

From an engineering standpoint, the ability to perform gate operations on any 
individual qubit or set of qubits with a given phase at any step in a series of operations 
is requisite to a universal quantum computer. We assume that all operations are 
synchronized to a local oscillator in perfect resonance with the qubits. Each qubit 
initially has an arbitrarily defined phase, and subsequent phase accumulations from 
interactions must be tracked so that each operation account for the phase of the 
individual qubit. These interactions are primarily AC Stark shifts from the optical 
control fields and Zeeman shifts from ambient magnetic fields. Our goal here is to 
prescribe a set of gates that leave the qubits independent of the optical phase of the 
driving field after each operation is complete while enjoying passive isolation from Stark 
and Zeeman qubit phase shifts. 

In this paper, we consider several quantum gate schemes in the trapped ion system, 
concentrating on a class of currently-favored quantum gates that rely on a "spin- 
dependent force" [HI EH EEl HU EH H2 HI] , and are relatively insensitive to motional 
heating from noisy background electric fields [TS1 Q2j- We will discuss the sources of 
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phase decoherence for these gates and describe methods to suppress decoherence from 
slow phase drifts of the driving optical fields. We find that certain gate schemes can be 
simultaneously insensitive to background magnetic fields, making them quite robust for 
long-term computations. 

This paper is divided into three sections: Section 1 lays out the background for 
the various quantum gate schemes for trapped ion system, covering detailed steps for 
performing single qubit operations and a discussion of the original Cirac-Zoller model 
with special attention paid to the phases of the qubits [Zj . Section 2 describes entangling 
gates using a spin-dependent force. We show how special arrangements of the classical 
driving fields can suppress slow phase decoherence from laser noise and external magnetic 
field noise j!7j . Section 3 shows how to similarly suppress long-term phase noise in a 
recent "fast" gate scheme that relies on impulsive spin-dependent optical forces |12j . 

1. Background 

Typical atomic ion species for quantum information applications such as 9 Be + , 43 Ca + , 
and 111 Cd + have a single valence electron with a 2 Si/2 ground state and 2 Pj first excited 
electronic states. In isotopes with nonzero nuclear spin, the 2 5i/2 ground states are split 
by the hyperfine interaction. An applied static magnetic field B provides a quantization 
axis and removes degeneracy in the ground state Zeeman levels. Two states, one from 
each ground state hyperfine level, are designated as the qubit states, denoted by fj and 
li for each ion i and separated by energy — E± = hw . At certain values of B , this 
energy splitting can be insensitive to magnetic field fluctuations to first order, forming a 
qubit that can have particularly good phase stability. Such qubit levels are termed "clock 
states" because their stability is exploited in atomic clocks [20]. The qubit frequency 
splitting ujq is usually in the microwave range and large compared with the radiative 
line-width 7 of excited electronic states, therefore the qubit can be measured with high 
fidelity by resonant pumping to a cycling transition between one hyperfine state and an 
excited electronic state jSJ. Initialization of the qubits can also be accomplished with 
high accuracy using optical pumping techniques (Figure EJ. 

We assume ions are confined in a linear Paul trap [21] with a combination of 
static and radio frequency (rf) electric quadrupole potentials [22]- When the ions are 
sufficiently cold, their Coulomb repulsion balances the external confinement forces, and 
the ions form stationary crystals with the residual motion described by coupled harmonic 
oscillatory motion in the trap. The number of collective modes scales linearly with the 
number of atoms N in the trap, making it difficult to isolate and control all modes of 
oscillation for large numbers of ions. To circumvent this difficulty, an architecture has 
been proposed for shuttling ions between multiple trap regions in a trap structure such 
that it is only necessary to localize a small number of ions at a given time 0. Two- 
qubit quantum gates in addition to arbitrary single qubit rotations are sufficient for 
the engineering of arbitrarily complex entangled states P , so we focus on the operation 
of quantum gates on N = 2 ions, although extensions to larger numbers of ions is 
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straightforward. 

The ions arrange themselves along the weakest (x) axis of the confinement potential, 
and the position operator of each ion can be written as Rj = Ro,j with the operator 
f j describing the small quantum harmonic oscillations of each ion about its equilibrium 
position Ro,i- Of the six normal modes of oscillation for the two ions, only the two axial 
normal modes will be considered for simplicity (see section 1.1): the center-of-mass 
coordinate q\ = (x\ + x 2 )/V^ and a "stretch" coordinate q 2 = (xi — x 2 )/\/2, where Xi is 
the component of f j along the a;-axis. The base Hamiltonian for the collective system is 

^0=5^^° 11*) (T*l + (i) 

1=1,2 v=l,2 

where u> is the frequency difference between the two qubit states; iQ\ and tu 2 = V3ui are 
the frequencies associated with the center-of-mass and stretch modes, respectively; and 
aj, and a v are their respective harmonic oscillator creation and annihilation operators. 
The first term in the Hamiltonian describes the internal energy of the ions, and the 
second term describes the external vibrational energy of the system. 

Single qubit rotations between hyperfine qubit levels (not involving ion motion) 
can be performed by applying appropriate radiation fields. For example, a resonant 
microwave field can directly couple the qubit levels through a magnetic dipole 
interaction, resulting in coherent Rabi oscillations between the two qubit states. 
Alternatively, optical stimulated- Raman transitions can be employed using two 
optical sources that coherently couple the qubit states through excited 2 Pj electronic 
states, as discussed next. 

1.1. Coherent Interaction between Trapped-Ion Hyperfine Qubits and Optical Fields 

An optical coupling between the hyperfine qubit states and an excited electronic state 
of each ion can be exploited to entangle qubit states with collective motional states, 
forming the backbone of most trapped ion quantum logic gates. As shown in Figure 121 
we assume each of two trapped ions consists of three levels: the two ground state qubit 
levels and and an excited electronic state |ej) having respective optical frequency 
spans u^e and o)j >e = Q^ e +uo. Two optical fields E;(r) = Ei(r) cos(k/-r— uit — 0/)e/ with 
I = a,j3 and polarization e h connect each of the qubit levels and respectively, 
to state |ej) through electric dipole operators and fi<. We assume that the optical 
fields have a difference frequency ojp — u a = u Q + 5oj and are both detuned from the 
excited-state resonance by A = oj^ e — u a , as drawn in Figure El These fields evaluated 
at the ion's position E^(Rj) are what ultimately couples the spin to the motion. 

The interaction can be transformed to a rotating frame at frequency u a in order to 
remove all terms varying with optical frequencies, and under the usual optical rotating 
wave approximation (RWA), the interaction Hamiltonian between the fields and the ions 
is 
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^ = \ X) [ - (^.a.ie*--* 4 "^ |e<) (1,1 + /i.e.) - (^^e^'^-^e" 1 ^ |e<) <T*1 + &-c) 

i=l,2 

+ A|e i )(e i |l (2) 



In this expression, the dipole coupling strengths between qubit state \rrii) = , \U) 
and excited state | e^) from laser field I on ion i are given by hg m ^^ = — \i m ■ e;^(R i )/2. 

For most of the remainder of this paper (outside of Section 3), we assume that the 
detuning A of the optical fields from electronic resonance is much larger than the excited 
state line- width 7 and the couplings \g m ,i,i\ , so that spontaneous emission during the 
optical coupling is negligible [S] and the excited state |ej) can be adiabatically eliminated. 
Applying RWA on the microwave frequencies, we find 



k ' = ~\ X [ (n ie -^-^^ \w (U\ + h. c ) 



h 

i=l,2 

+ xu\U)(U + x^\h)(h\ (3) 



where Ak = kg — k a and A0 = <fip — <p a are the differences in the wave-vector and the 
phase of the two applied fields, Qi = g^ >a) ig : f /3i /2A is the "base Rabi frequency" directly 
coupling the qubit states of ion i, and Xm,i = (\g m> a,i\ + \9m,p,i\ ) /2 corresponds to the 
AC Stark shifts of the qubit level \rrii) of ion % by both optical fields. 

For simplicity, we assume Ak is parallel to the x-axis (Ak = |Ak|), so that the 
interaction deals only with axial motion (although it is straightforward to treat the 
more general case). We substitute the x-component of the position operator Rj for ion 
% by 

X, = X 0ti + ^={a x + a\) ± ^={a 2 + a\) (4) 

where q v = a/ hj (2Mu u ) is the root mean square spatial spread of the ground state 
wavepacket for normal mode v of oscillation in the trap, M is the single ion mass, 
and the plus (minus) sign refers to ion % = 1 (i — 2). In the interaction frame of the 
vibrational levels, Equation El becomes 

^ i=l,2 

+ x/3,i\U){U\+x a ,i\W(1i\] (5) 

The Lamb-Dicke parameters are defined by 771 = Akqi/ \pl and r] 2 = Akq 2 / v2 = r\\j v3j 
representing the coupling strength between the fields and each normal mode. 

The above treatment can be generalized to the case of multiple optical sources 
that connect both qubit states to any number of excited states, resulting in higher 
order expressions for and Xm,%- Here however, we are mainly interested in the 
sensitivity of entangling gate operations on the optical phases 4>i. The net optical 
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phase appearing in the coupling Hamiltonian (Equation |HJ) is sensitive only to the phase 
difference A<ft = <pp — <p a between the two optical fields, so that when both fields are 
generated from a single laser and modulator, fluctuations in the optical phase of the laser 
source become common mode and do not lead to decoherence [21]. However, in order to 
couple the qubits with the motion for entangling quantum gates, the optical sources are 
generally non-copropagating (Aky^ 0), opening up the sensitivity to decoherence from 
fluctuations in relative beam path lengths or ions' positions through the phase factor 
e iAkx 0ii -iA(/> ^ This requires interferometric stability between the optical paths of the 
fields E a and Ep, which should be feasible over short times using stable optical mounts 
and indexing the laser beams to the trap structure itself. However, over the long time 
scale represented by an extended quantum computation, drifts in the phase e lAkX o,i-^<t> 
can be a serious source of decoherence. 

1.2. Resolved-Sideband Limit 

Equation |3] includes direct couplings between qubit states, and entangling couplings 
between qubit states and trapped ion motional states. We consider the case where 
the base Rabi frequencies f2j are much smaller than the vibrational frequencies uj u of 
the ions in the trap. In this case, the difference frequency 5u> of the optical sources 
can be tuned to particular values so that a single stationary term emerges from the 
above Hamiltonian, and all others couplings can be neglected under the rotating wave 
approximation. In this regime, as seen from the rest frame of the ions, the applied 
laser fields acquire resolved frequency- modulation sidebands from the ions' harmonic 
vibration. We concentrate on three spectral features: the "carrier" , the first upper and 
lower sidebands, each selected by appropriate tuning of the radiation field difference 
frequency. 

1.2.1. The carrier When the difference frequency between the optical sources is tuned 
to the free-ion qubit resonance (compensating for possible differential Stark shifts, 
assumed to be equal for the two ions), then uip — u) a = u + xi,i — Xi,i = (Figure 
) , and we find the stationary term in Equation El is given by [H] 

Hi ar = ~ EE {^V nun2 e^ kX ^~^af + h.c) \n u n 2 ) (n u n 2 \ . (6) 

i=l,2 ni,ri2 

This "carrier" interaction describes simple Rabi flopping between the qubit states in each 
ion, where the qubit raising and lowering operators are defined by o+ = (U\ and 
= ( ^))t = ||.) ( T .|. Also, V nun2 = e-^lW 2 )c ni ^l) Cn2 ^l) is the Debye-Waller 
factor that exponentially suppresses the carrier coupling due to ion motion described 
by vibrational quantum numbers n v in each mode v of motion, with C nv (z) being a 
Laguerre polynomial of order n v When the ions are confined to the Lamb-Dicke 
limit (LDL) where rj%(n u + 1/2) 1, then T> nijTl2 ~ 1. However, if the motion of either 
mode is not in a pure eigenstate of harmonic motion and outside of the LDL, then the 
Rabi frequency will depend upon the noisy motional quantum state and lead to qubit 
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decoherence. In order to avoid this problem, carrier operations are often performed with 
copropagating Raman beams (j] u = 0) thereby forcing T> nitU2 = 1. A carrier transition 
using a copropagating geometry can also be insensitive to the phase noise of the source 
laser, making it ideal for single qubit rotations. 

1.2.2. The first lower (red) sideband When the difference frequency between the optical 
sources is tuned lower than the free-atom qubit resonance by the vibrational frequency 
uj v of mode v (again compensating for differential Stark shifts), up — uj a = oj' — u„ 
(Figure Eb) and we find the stationary term in Equation El is |Hj 

n? b = ~E (^e i(AfeXo '*~ A ^U^ + h - c ) ■ ( 7 ) 

i=l,2 

This "red sideband" interaction describes Rabi flopping between the coupled 
qubit-motional states \{,n u ) and — 1) in each ion, where T)' nvn t = 

e-Uvi+v*) Cn ^~ 1 ^ £>n v , (vl') i s the Debye- Waller factor for the first sideband, with v' ^ u 
the "spectator" mode of motion. Here, £*(z) is an associated Laguerre polynomial. 
This interaction is the Jaynes-Cummings Hamiltonian j2S] where energy is exchanged 
between the internal qubit and the external harmonic oscillator states. 

1.2.3. The first upper (blue) sideband When the difference frequency between the 
optical sources is tuned higher than the free-atom qubit resonance by the vibrational 
frequency u u of mode v (once again compensating for differential Stark shifts), ujp— u Q = 
uj f Q + oj v (Figure Efc) and we find the stationary term in Equation El is now jHj 

^ = ~\ E (^^^-^)^ n ^>at + h.c) . (8) 

i=l,2 

This "blue sideband" or anti- Jaynes-Cummings interaction describes Rabi flopping 
between the coupled qubit-motional states \i,n v — 1) and \],n u ) in each ion. 

1.3. The Cirac-Zoller Gate 

The original Cirac-Zoller (CZ) scheme [TJ I2H1 illustrates how entanglement between 
trapped ion qubits can be achieved through coupling of each qubit to a common mode 
of motion in the trap. The CZ scheme allows the operation of a controlled-NOT gate 
between two trapped ion qubits, flipping the state of a target qubit (e.g., || 2 ) ^ IT2)) 
only when the control qubit is, say, in state This can be accomplished by cooling 
a collective motional mode v of the two ions to the [0^) ground state and performing 
the following three steps: 

(1) A carrier 7r/2-pulse on the target qubit with associated phase 0, 

(2) A 7r phase gate on two ions 

(3) A carrier — 7r/2-pulse on the target qubit with phase (step (i) reversed). 
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ITT) - 


► ITT) 


iu> - 


> IU> 


lit) - 


> III) 


III) - 


> -III) 



Steps (1) and (3) are simply carrier couplings on the target qubit ion, achieved by 
focusing radiation on the target ion only (f^ = 0) and applying the radiation for a time 
tw/2 (|^2^7r/2 = 7r /2). Step (1) results in the evolution 

« IT2) + /3 1| 2 ) - 1 ^ ^ IT2) + ^ ^ ^ II2) • (9) 

Step (3) is identical to step (1) except the phase is shifted by ir. There are many ways to 
implement step (2), one of which is: i) a 7r-pulse blue sideband that maps the internal 
qubit state of the control qubit to the collective state of the ion pair, ii) a 27r-pulse 
coupling the || 2 ) \n = 1) state exclusively to an auxiliary level, and iii) a 7r-pulse on the 
blue sideband to map the collective motional state back to the control bit (see Figure 
HJ). The net effect of these steps produces the following phase gate: 



(10) 



Here every state maintains a constant amplitude and the phase is well-defined. However, 
steps (1) and (3) contribute an additional phase to the controlled-NOT gate: 

ITT) - ITT) 

ITI) - |U> (u) 
lit) - e*|U> ' 1 ] 

III) - e^UT) 

The phase of the Cirac-Zoller controlled-NOT gate therefore depends solely on the phase 
of the 7r/2-pulse single qubit rotations. As mentioned in section 1.1, the sensitivity of 
single qubit rotations to optical phase can be removed using copropagating Raman 
beams requiring only a stable microwave source driving an optical modulator. This 
conversion between a phase gate and a controlled-NOT gate is extremely useful for 
many entangling gate schemes, since phase gates are more intuitive to construct and 
have an inherently well-defined phase. 

The CZ model for trapped ion quantum logic gates has many drawbacks, including 
the need for individually addressing the ions with optical sources, and the requirement 
that the motion be prepared in a pure state of collective motion, usually through laser- 
cooling to the |0„) state. In the remainder of this paper, we consider improved schemes 
for trapped ion quantum gates that do not have these requirements. In some cases, we 
will see that the sensitivity to the optical phase e lAkX o,i-i-^ can a i so b e suppressed. 



2. Spin-Dependent Forces in the Resolved-Sideband Limit 

Unlike the Cirac-Zoller gate where the internal qubit state of one ion is directly 
transferred to particular eigenstates of motion, entangling gates using spin-dependent 
forces coherently displace the initial motional state in the position/momentum phase 



Phase Control of Trapped Ion Quantum Gates 



9 



space, a process through which each spin state can acquire an independent geometric 
phase [HI EU Ell EH! S3 ■ The nonlinearity in these phases can result in a final state that 
can no longer be separated into two independent qubit subspaces, thus entangling the 
internal states of the two ions. This produces a phase gate similar to the Cirac-Zoller 
scheme, which can then be converted to a controlled-NOT gate when combined with 
single qubit rotations. 

In this section, we focus on gates in the resolved-sideband limit where the interaction 
time is much longer than the trap period. The interaction Hamiltonian is proportional 
to &i ■ ni <g) & 2 ' n 2> where <x, is the Pauli spin matrix operating on the internal qubit 
states, and rij is a unit vector pointing in a particular direction on the Bloch sphere 
for ion i. The eigenstates of &\ ■ n x £g> cr 2 ■ n 2 each experiences a different force from 
the interaction (see Figure EJ). The gates are categorized according to the direction of 
rij: in a "<r 2 gate", the differential force is applied via a differential AC Stark shift on 
the states and ||,) induced by the laser fields [T£]. However, clock states exhibit no 
differential AC Stark shift when the Raman detuning A is large compared to the qubit 
frequency splitting u (see Appendix A) , so the only available qubit states for a b z gate 
are thus susceptible to magnetic field fluctuations. In a "o^ gate" , optical fields driving 
spin flips and coupling to the motion produces a differential force between eigenstates 
of &i ■ (p^ where the unit vector (p i = cos(0j)x + sin(0j)y lies on the equatorial plane of 
the Bloch sphere [T H HT H l2T | ITT] . Although this gate is compatible with clock states, the 
optical beam configuration can give rise to extreme sensitivity of the qubit phase on the 
optical phase of the driving field, which can be the limiting factor in the fidelity of the 
gatejUJ. In this section, we propose a method to cancel this phase dependence on the 
optical field, relaxing the constraint on long-term interferometric stability between the 
two Raman beam paths for the entirety of a multi-gate sequence quantum algorithm. 



2.1. Forced Quantum Harmonic Oscillator 

In order to understand the spin-dependent force, we start by considering the effects 
when a force is applied to a harmonic oscillator. In general, a forced harmonic oscillator 
has a Hamiltonian of the form [2E] 

H = hu(a)a + -) + f*(t)x a + f(t)x a\ (12) 

where a and a) are the annihilation and creation operators respectively, and Xq = 
y/KJ (2Mw) is the root mean square spatial spread of the ground state wavepacket. The 
first term is the unperturbed Hamiltonian for the harmonic oscillator of frequency u, 
and the last two terms correspond to an external time-dependent force f(t) applied to 
the system. In the interaction picture 

H^t) = f*(t)x ae- iut + f(t)x Q tfe itJt . (13) 

Assuming the force f(t) = Fe~ l ^~ s ^ /2 is detuned from resonance by frequency 5 to, 
then the interaction Hamiltonian can be rewritten as 

Hl (t) = ^ae-™ + ^aV« (14) 
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The state after an interaction time t is prescribed by the time-evolution operator 

U(t) = exp | -I ^Hj(t')dt' + ±J*dt'J* dfiHjit'),^")} + ...J | .(15) 

If we consider only the first term in the exponent of the evolution operator and 
substituting in the interaction Hamiltonian from Equation El the resulting operator is 
exactly the displacement operator 

D(a) = e aat+Q * a , (16) 

with a defined as 

1 f tFx 0A8t> 



a(t) = -- J -^e M dt'. (17) 

The displacement operator translates motional states in position/momentum phase 
space without distortion (Figure EJ). For example, a displacement on an initial ground 
state of motion results in a coherent state |a) = D(ct) |0), where the final state is defined 
in terms of number states as \a) = e~^\ a \ Yl^Lo ^7=? \ n ) ■ I n terms of x-p coordinates, 
a = (l/2x )(x + ip/Mu). 

The remaining higher order terms in the time-evolution operator originate from 
the non-commutative property of the interaction Hamiltonian at a given time with 
itself at different times. This can be understood by considering the displacement 
operators, which do not commute with one another but rather follow the commutation 
rule D(a)D(fl) = D(a + (3)e %Im ^ a,3 *\ Therefore the complete time-evolution operator 
can be constructed by integrating over infinitesimal displacements in time: 

U(t) = e iq,{t) t)(a{t)), (18) 

with the geometric phase accumulated over the entire path from time to t expressed 
as 

$(t) = Im{ [ a{t')*da{t')). (19) 
Jo 

For a near-resonant driving force with detuning 5 ( Equation ITlj). the initial motional 
state moves in a circular trajectory of radius F/(2hS) with periodicity T = 2ir/5 in the 
rotating frame of harmonic motion, following the path (from Equation ITTj) 

In one period of evolution under this force, the motional state returns to its original 
phase space coordinates, but acquires a geometric phase of 

% = 7T|FX ° |2 (21) 
2(M) 2 K ' 

equivalent to the area enclosed by the trajectory (Figure EJ). 

For a single qubit experiencing a spin-dependent force, the interaction Hamiltonian 
includes a dependence on the internal spin state of the ion. Assuming the force 
couples to only one of the vibrational modes and the other mode can be neglected 
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under the rotating wave approximation, the most general expression for the interaction 
Hamiltonian is 

Hi= (^pae- iat + ^paV a A |m> <m| , (22) 

where m denotes the internal qubit state of the ion, and || n ) and || n ) are the eigenstates 
of a • n associated with eigenvalues +1 and —1 respectively. The interaction provides 
no coupling between the two orthogonal spin eigenstates of <x • n, but the motional state 
becomes entangled with the spin state as the differential force pushes the motional states 
of the two spin components in separate directions. At time t = 2nir/5, where n is an 
integer, the two motional states overlap again, disentangling the vibrational component 
of the wavefunction from the spin, but leaving the spin component with a phase shift 
due to the difference in the geometric phases of the paths. The interaction Hamiltonian 
can also be written in terms of the <x • n operator as follows: 

Hj = f^ae- i5t + *jpat e «*) J+(^pae-^ + ^tfjA <xn, (23 

where / is the identity operator, F + = (F^ n + F^ n )/2 and F + = (F^ a — Fj n )/2. 

Now consider a spin- dependent force applied simultaneously to two ions in the same 
trapping potential. The total force on the system is now dependent on the spins of both 
ions. The interaction Hamiltonian now becomes 

H I= J2 ( Fml ' m f )X ° a^ + F ^ t)XO a\ \ mhm 2) ( m l, m 2\ , (24) 

ml,m2=Tn,J.n V " -J 

where mi, m<i denotes the internal qubit state of ion 1 and ion 2 respectively, and 
F m i,m2 = Fmi + F m 2 is the total force applied to the state |ml,m2). The geometric 
phase of an enclosed loop is proportional to \F mitm2 xo\ 2 /5 (Equation EH), which can 
be calibrated so that the nonlinearity results in a wavefunction whose spins are not 
factorizable, thus creating entanglement between two ions. 

The following sections will provide specific examples of entangling gates using spin- 
dependent forces. While the fundamental concept is the same in both instances, the 
experimental requirements and the susceptibility to various sources of phase decoherence 
are distinct. We will discuss these cases in detail and provide some solutions for phase 
control of these gates. 

2.2. The a z - gate 

As the name implies, the <r z -gate applies a differential force on the eigenstates of the 
unperturbed Hamiltonian. The interaction has no coupling between the two eigenstates 
of <t 2 , conveniently avoiding the neccessity to produce large frequency shifts in the 
Raman beams to bridge the hyperfine splitting. Instead, the Raman beams couple to 
the vibrational states without driving a spin flip. The two Raman beams form a beating 
wave with periodicity 2irc/ {oj v — 5) along the weakest trap dimension. Due to the AC 
Stark effect, this wave form a moving periodic potential, exerting a near- resonant force 
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on the ion in the direction of travel. If the AC Stark effect has different amplitude on 
the two qubit states, then the two states experience different forces from the beating 
wave [IB]. 

The a z gate uses two non-copropagating beams with frequency difference u v — 5, 
where uj v is the frequency of vibration and 5 is the detuning from the vibration frequency. 
For this example, we will let the beams couple to the stretch mode u = uj 2 = y/3ui, 
though the same algebra can be carried out for the center-of-mass mode. (Stretch 
mode is a better candidate since it exhibits lower levels of decoherence from background 
electric fields |IB|.) We apply two fields S A e < ^ ,x - WA *-*^eA+SBe < ^ fl,x ~ UB *~* B) eB where 
the frequency difference ujb — oja = ^2 — 8 is slightly detuned from the stretch mode 
frequency. The field couples each of the spin states to the excited P state, and is detuned 
by a large frequency A (see Figure HJ). Using the same RWA and adiabatic elimination 
of the excited state used to obtain Equation El the interaction Hamiltonian for a single 
ion becomes 

Hi = ^ { [M 2 + l^l 2 + (2 9 ; A9 , B e^-^- s ^ + h.c)] | T > (T| 



+ 



\9lA\ 2 + \9 lB \ 2 + (2gl A g lBe ^-^- s ^ + h.c.)] ||> (||} (25) 



where g m> i = /J, m ■ eiEi/2h is the single photon Rabi frequencies associated with each 
field I coupling qubit state \m) to excited level |e), Ak = k B — k A is the wave vector 
difference, and A0 = 0# — <Pa is the phase difference between the driving fields. The 
first two terms in the expressions in Equation EH1 contribute to the average Stark shift for 
each of the states, which can be canceled by carefully choosing the polarizations e x and 
e 2 The cross terms represent the variation in intensity formed by the interference 
pattern that pushes the ion, and must have a different magnitude and/or phase between 
the two qubit states to create a differential force. In the Lamb-Dicke limit, and assume 
the detuning A is approximately the same for both spin states (A 3> a> ), the interaction 
Hamiltonian for two ions can be written as 

i=l,2 mi =1,l 

( — Y^ a ^ e + — 4 e J \mim 2 ) {m 1 m 2 \ , (26) 

mi,m2=T,l 

where F mum2 q 2 = (hrj 2 D' n2 n J 'A)(g^ l!A g mi>B e i<t ' 1 -g^ 2tA g m2; Be i ' t>2 ), and ^ = AkX , t -A(j). 

The phase difference between the force applied to the two ions is determined 
by the optical phase difference <pi — <p 2 , which corresponds to the ion spacing at 
equilibrium. If the ions are spaced by an integer multiple of the optical wavelength, 
i.e. Afc(Ao i i — Xq, 2 ) = 2nn, then they experience the same phase in the force, i.e. 
0i — 02 (see Figure |HJ). This is a convenient case since the forces cancel when the two 
spins are aligned in the same direction, and displacement occurs only when the spins 
are anti-aligned. A physical explanation of this scenario is that the stretch mode can 
be excited only when the two ions are pushed in different directions or with different 
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magnitudes. The fastest gate time possible for this scheme is when the anti-aligned 
states acquire a tt/2 phase shift in time T = 2tt/5, or in other words, a round-trip 
geometric phase $o = tt |F|^g 2 | 2 /2{h5) 2 = ir/2. Under these conditions, and assuming 
all the average Stark shifts have been accounted for, the gate performs the operation 

PUD!!: 



ITT) - 


- ITT) 


in) - 


- i|U> 


ur> - 




in) - 


- in) 



(27) 



With a phase shift of —tt/2 on both qubits, the final state is equivalent to the result 
from a standard phase gate in Equation HH 

Note that the end result is completely independent of the optical phase of the drive 
field. The optical phase fa is absorbed in the term F mi>m2 , translating to a phase shift in 
a that defines the coherent state. Since the acquired geometric phase depends only on 
the area enclosed by the trajectory, the phase of the resulting state has no correlation 
to the phase of F mitm2 . While the optical phase still needs to be coherent during a 
gate, variations in phase between gates are acceptable since they have no impact on the 
outcome. 

Ideally, if the a z gate can be performed on magnetic field insensitive states, the 
gate would be extremely phase stable. Unfortunately, in the limit where the detuning 
from the excited state is large compared to the hyperfine splitting, magnetic field 
insensitive states have no differential Stark shift (see Appendix A). Therefore, this 
gate is susceptible to decoherence from magnetic field fluctuations if performed on 
magnetic field sensitive states, or if performed on magnetic field insensitive states in 
the limit where detuning from excited state is comparable to the hyperfine splitting, by 
spontaneous emission from the excited state. 

2.3. a 0- gate 

We call the gate scheme proposed by Maimer and S0rensen (HUH] a ^4> S a t e because the 
interaction is analogous to the a z gate operating in a rotated basis. Although the original 
treatment describes the interaction in a four-level ladder system, here we describe it in 
terms of spins and displaced motional states as in section 2.1. The M0lmer-S0rensen 
gate employs simultaneous addressing of both ions with bichromatic fields, one detuned 
from the red sideband of a vibrational mode by frequency 5 and the other from the 
blue sideband by —5. The two sidebands have equal strength i] u Q/2 in the Lamb-Dicke 
Limit, and once again we assume the force couples only to the stretch mode. The 
interaction Hamiltonian is simply the sum of the red sideband plus the blue sideband 
(from Equations [7| and (HJ) with a detuning 5: 

Hi = -^J2 ^fi l ^ 2 ,^(e j(Afc " Xo -- A *" ) ^ ) a 2 e- i<5 * + e ^hx 0li -A M & yst + h c ^ (2g) 

i=l,2 
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where r] 2 Q Q D' n2 n , is the sideband Rabi frequency, Ak r and A&& are the wave-vector 
difference for the red and blue sidebands, X 0) i indicates the equilibrium position of the 
i-th ion along the x-axis, and A(j) r and A<pb are the phases of the red and blue sidebands 
respectively. We can simplify this expression to: 

i=l,2 

2^ ( — -j^a 2 e lM + — if^-ale J \mim 2 ) {m x m 2 \ , 

(29) 

where 

<7« . = <Tj • [cos(<^)x + sin(0 5 ,i)y] = <^V^ + o^e***. (30) 

Here Fi = Hr] 2 Q,iD' n2 n ,J q 2 is the differential force on the z-th ion, 4>s,i — — (Ak r X 0t i — 
A(p r + AkbX 0ji — A(p b )/2 is the spin phase of the i-th ion, (j> M>i = (Ak r X 0ii — 
A(p r — AkbX 0ji + A(p b )/2 is the phase of the force on the i-th ion, and F mijma = 
±Fie~' u ^ M ' 1 ± F 2 e~^ M ' 2 where ±Fj corresponds to the force on the spin state =J \^, Si 
, X^s- i respectively on the i-th ion. As in the & z gate, we set the phase of the force acting 
on the two ions to be opposite, i.e. Fie % '^ M ' 1 = —F 2 e^ M < 2 , and choose 5 and F such that 
the round-trip geometric phase satisfies the condition $ = 2tt \Fiq 2 \ 2 / (hS) 2 = n/2. 
Then the final state of the gate is equivalent to the final state in Equation EH except 
|fi) and the eigenstates of &z , are replaced by |t<fe,i) anc ^ |-Us,i)> the eigenstates 
of . . This gate written in the a z basis produces the following truth table: 

ITT) - ^{ITT)-^-+^)|U>} 

III) - ^{|U>-iUT» r . 

UT> - ^{|it>-i|U» 1 j 

III) - ^{IU)-ie-^+^)|TT)} 

Note that after the gate, only the spin phase remains, while the motion phase has no 
effect on the final state. As in the a z gate, drifts in the motion phase between gates is 
acceptable. However, the spin phase must be maintained between gates, or alternatively, 
an equivalent entangling gate with the dependence on the spin phase removed can be 
formed using a combination of gate and other quantum operations. 

An analysis of noise sources for the spin phase requires careful consideration of 
the physical experimental setup in the laboratory. To drive the red sideband and 
the blue sideband transitions simultaneously, a minimum of three optical frequencies 
are required, assuming one frequency can be used for both sideband couplings (see 
Figure El). Since each of the two fields driving a sideband must have a non-zero wave- 
vector difference Ak, the optical beams can be set up such that two fields at different 
frequencies share the same wave- vector kg, and each individual field combined with 
a third field having a different wave-vector k^ drives a sideband transition. In other 
words, if the field along k^ has frequency lja, then the field propagating along kg needs 
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both a frequency component uja ± (uj' — uj 2 — 5) to drive a detuned red sideband, and a 
frequency component uja ± (uj' + u 2 + 5) to drive a detuned blue sideband. The choice 
of the positive or negative frequency differences between the fields determines the sign 
of Ak r and Ak b , and determines the gate's susceptibility to the phase stability between 
the two wave-vectors. 

2.3.1. Phase sensitive geometry The first scenario involves frequencies of both fields 
along k B being higher (or lower) than u>a, then the wave- vector difference Ak for both 
the red and the blue sideband propagates in the same direction (Figure ITUk). For 
example, let the field along k B include both u>a + u>' Q — lo 2 — 8 and u>a + u' + u> 2 + 8 
frequency components. Then the wave- vector difference Ak r = kg — = Ak& for 
the red sideband and the blue sideband point in the same direction. Instability in 
the relative beam paths results in an equal phase shift in the sideband transitions, i.e. 
5(f> r = 5(f>b = S(f>. This results in a net shift in the spin phase by 84>s,i = S(j>- This is not 
a desirable situation since the outcome of the gate is sensitive to changes in the beam 
path length difference on the scale of an optical wavelength. 

However, we note that the spin phase shift is exactly the same as the phase shift on 
the non-copropagating carrier transition, i.e. when the field propagating along k B has 
frequency u>a + ^o- Therefore it is possible to construct a phase gate using the following 
Ramsey experiment: 1) Perform a ir/2 rotation on both ion with phase shift 5<ps,i = 
using the non-copropagating transition; 2) Perform the 0$ gate using the frequencies 
listed above; 3) Perform a —tt/2 rotation on both ions with phase shift 5<fts,i = $4> using 
the non-copropagating transition. This rotation from z to before the gate and the 
subsequent rotation back to z after the gate effectively removes the dependence on 
the spin phase <fi as long as the spin phase is constant during the Ramsey experiment. 
The final state becomes identical to Equation |2H and has no residual dependence on S(j). 
In addition, this scheme is also insensitive to ion spacing since the phase of the push 
force is always zero in the basis defined by <ps,i- 

2.3.2. Phase insensitive geometry Another scenario is to select the frequencies along 
k B to straddle the frequency along k^. Then the wave- vector difference for the red and 
the blue sideband propagates in opposite directions. For example, let the field along k B 
include both oja — — uj 2 — $) and uja + uj' + uj 2 + 5 frequency components (Figure 
[TUb). Then the wave- vector difference for the red sideband Ak r = — kg + k^ = — Ak& 
is in the opposite direction as the blue sideband. Instability in the relative beam paths 
results in an opposite phase shift in the sidebands, i.e. — S(j) r = S(f>b = Sep. This results 
in a net zero change in the spin phase 5<ps,i — 0, removing any spin phase dependence 
on 5(f) from the gate. Hence this configuration is termed "phase insensitive" . 

However, the motion phase in this setup acquires a dependence on the phase shift 
5<p. Therefore the phase of the force on each ion should be calibrated to be the same 
by setting the ion spacing (using the trap frequency as a tuning parameter) equal to 
X 0j i — X 0j2 = 2mr/Ak, where n in an integer. While it is possible to produce similar 
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entanglement operations with other values of ion spacing, the gate speed will be slower 
for the same intensity from the laser, and the output will appear different than the 
expression in Equation 13 II 

Similar to the other configuration, it is possible to construct a phase gate with the 
transformation in Equation |2H using an analogous Ramsey experiment involving single 
qubit rotations in phase with the gate: 1) Perform a it/ 2 rotation on both ion with 
phase shift 5<ps,i = using either a calibrated and phase locked microwave source or 
a copropagating Raman transition; 2) Perform the gate using the frequencies listed 
here; 3) Perform a —it/2 rotation on both ions with phase shift 5<ps,i = 0. 

3. ovgate with fast pulses 

The <r z -gate can also be achieved by applying spin-dependent momentum kicks on the 
ions with fast laser pulses [121 EHj • For gates in the resolved-sideband limit discussed 
in section 2, the ion is assumed to be confined within the Lamb-Dicke limit, where the 
spread in the position of the ions from their equilibrium positions is much smaller 
than the optical wave length. Outside of this limit, the effective Rabi frequency 
fluctuations leads to significant gate errors. For gates using fast laser pulses [1211201, the 
impulsive spin-dependent force from the traveling wave has an almost uniform intensity 
distribution around the ion's position, and this fluctuation of the effective Rabi frequency 
can be safely neglected. These pulsed gates can therefore faithfully operate outside of 
the Lamb-Dicke limit. In this section, we want to show that non-trivial phase errors can 
arise when the ions are outside the Lamb-Dicke limit, and suggest a method to cancel 
these errors by carefully selecting the direction and timing of momentum transfer, a 
technique reminiscent of the phase cancellation effect in the phase insensitive gate 
configuration discussed in section 2.3.2. 

The central component of the fast a z gate in the context of ground state hyperfine 
qubits is a set of fast resonant laser pulse pairs that exclusively couple one of the two 
qubit states (here taken to be ||)) of each ion to the excited state |e). The coupling 
Hamiltonian follows from from Equation 2 with A = 0: 

Hj = hJ2- ^^e ifc - Xo ^-^e- ir?l(ai+a i )Tir?2(a2+a 2) | ef ) (|| + h.c^j (32) 

where g(t) is the resonant Rabi frequency of the transition for each ion and as before, the 
plus (minus) sign refers to ion i = 1 {i = 2). The pulse pairs are set to drive successive 
7r-pulses (Jq g(t)dt = %) from — > |e«) — > on the electronic transitions of each 
ion, with the pulse duration r taken to be much shorter than the radiative lifetime of 
|e) as well as the trap period 2ir/ui v . 

When these two successive fast pulses have non-copropagating wave-vectors and 
kg, and both 7r-pulses are completed in a time much shorter than the lifetime of state 
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\ei), the result is the following evolution for two ions 112j: 
T1T2 \ai) 2 -> T1T2 |a} 2 

T1I2 |a) 2 -> T1I2 |a + ^i)i|a-^2} 2 /ggx 
I1T2 |a) 2 -> I1T2 I a + %) 1 |a + ^ 2 ) 2 
Iil2|a)i|a) 2 I1I2 |a + 2%) 1 |a) a 

In this expression, |ai) and |a 2 ) are initial coherent states of the two modes of motion, 
and r\ v are the Lamb-Dicke parameters of the two modes associated with the wave- vector 
difference Akj = k^ — kg, exactly as defined in section 1.1. 

In the fast a z gate, a series of pulse pairs is applied to the ions so that the motional 
states of both modes of motion simultaneously return to the same position in phase 
space regardless of the state of the two qubits. When these fast pulses are interspersed 
with periods of free evolution of the two modes of harmonic motion, the result can be 
a a z phase gate for appropriate choices of pulse timing [121 EH] . This fast gate works 
independent of the motional state and outside of the Lamb-Dicke limit, as long as the 
motion remains harmonic. 

However, outside of the Lamb-Dicke regime, we find that this gate can be sensitive 
to changes in the phase of the optical fields due to the change in the position of the ions 
at different times. In order to see this effect, we note that this fast o^-gate involves a 
spin-dependent force on the ion from absorption of a photon from a laser pulse traveling 
in the k^ direction and emission of another photon to a pulse propagating in a different 
kg direction. We can lock the relative phase of these two propagating laser beams so that 
their phase difference is set to zero at the ion's equilibrium position. Then, if during the 
above impulsive kicks, the ion is at a position r from its equilibrium site, it will acquire 
a net spin-dependent phase factor of e lAk r from the absorption and the emission of the 
photon. This phase factor from each spin-dependent kick is non-negligible if the ion is 
outside of the Lamb-Dicke limit. For a complete gate operation with a series of laser 
kicks, with the spin-dependent phase factor for the j-th kick is denoted by e 1 rTj , the 
total phase factor after N kicks is given by e tipt with ip t = J2f=i (Akj • ij). If the gate 
speed is comparable or slower than the local ion oscillation frequency, the ion's position 
Tj at different laser kicks changes depending on the initial momentum and position are 
become almost uncorrelated. Therefore, the above effect contributes a random phase to 
the spin, which is a source of the gate infidelity. 

To eliminate this random phase effect when the ion is outside the Lamb-Dicke limit, 
one needs to require the gate speed to be significantly faster than the local ion oscillation 
frequency. In that case, the ion's position at different laser kicks are almost the same 
although they are still unknown. For any two positions rj and during the j-th and 
k-th kick respectively, the difference between them can be bounded as |r& — Tj\ < vT g , 
where v is the ion's typical speed and T g denotes the gate time. Due to this position 
correlation and the fact that the total momentum kicks Ylj=i = f° r the ^ as ^ g&te, 
we conclude that the random phase ip t is bounded by cp t < |Akj|t>T g , and the gate 
infidelity 5F = 1 — F from this random phase scales as (|Ak|t>T 9 ) 2 . The scaling of 
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8F can be further improved to (|Ak| vT g ) 2n if we use a more involved sequence of the 
kicking forces with n basic cycles. The gate time must be short enough to make the 
scaling parameter |Ak| vT g < 1. Under that condition, the gate infidelity can then be 
reduced rapidly to zero with an appropriate pulse sequence even if the ion is outside of 
the Lamb-Dicke limit |29j . 

Conclusion 

Most quantum logic gate schemes for trapped ions operate through interactions with 
optical electromagnetic fields. Some schemes, such as the gate and the fast a z gate, 
have a phase dependence on the phase of the optical driving field, which can become a 
major source of decoherence if uncontrolled. We have shown here methods to remove 
this phase dependence for these two entangling gates by choosing appropriate wave- 
vector orientations and pulse timings that naturally cancel the phase factor e lAk r upon 
the completion of the gate. Furthermore, the sideband resolved gate can operate on 
magnetic field insensitive qubit states, removing an unavoidable vulnerability of the a z 
gate. These techniques eliminates the random phase from the optical driving field while 
maintaining phase coherence at the rf or microwave atomic frequencies, allowing long 
gate sequences to be performed over the timescales beyond the coherence time of the 
optical fields. 
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Appendix A. Magnetic field insensitivity 

In this appendix we will show that magnetic field insensitive states have no differential 
Stark shift in the limit where the detuning from the excited state is much larger than the 
hyperfine splitting, i.e. A^p/A — > [31]. To find the field insensitive states for a system 
in the S1/2 ground state with some nuclear spin J, we write down the Hamiltonian for 
the hyperfine interaction in the presence of a magnetic field B: 

H = fi B ■ B + AI ■ J = gjB J + #/B I + AI J, (A.l) 

where J is the total angular momentum of the electron, I is the nuclear spin, and AI ■ J 
is the contact term, gj and gj are the Lande g-factors for the nucleus and the electron. 
The eigenstates of the Hamiltonian are linear combinations of the mp states, and can 
be represented as |^) = aAg] mj = ^,mj = mp :i — \g; mj = — |, mi = mp :i + |). 
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The coefficients a and b are functions of the magnetic field. If two states j^i) and j^) 
are magnetic field insensitive, then 



0. 



Applying Ehrenfest's theorem 



95 



r</* + gih) 



9 J , / ^ 



+ l&i 



(A.2) 



(A.3) 



Normalization of the eigenstates and solving Equation IA.2I gives the result |ai| = 
| a.2 1 2 H- giAuiF / '{gj — gi)- Since the dipole moment of the electron dominates the dipole 

« 10 -3 , 



moment of the nucleus, i.e. gi/gj 

|2 



1 |2 
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]a 2 | 


\U I 2 


N 



we can approximate it as 



(A.4) 

Now consider the Stark shift for each of these magnetic field insensitive states. The 
AC Stark shift is given by 

E ■ d\e, mjmj) (e, mjmi\ E ■ d 



Xi 



E 



mj , mi 



A - E 1 + Ei 



(A.5) 



where \e; mj, mi) is the excited state with the corresponding z-component of the electron 
and nuclear spins. Since the electric dipole only couples the orbital angular momentum 
of the electron, only couples to the states with the same mi. So the expression can 
be simplified to 



Xi 



E 



(g;mj = || E ■ d \e,mj)\ 
A-Ei + Ei 



+ N 2 E 



[g\mj = -|| E ■ d \e,mj)\ 
A-E 1 + E i 



, (A.6) 



mj mj 

If the energy difference between the two states ^1 and ^2 is small compared to A, and 
applying the results from equation IA.4[ then we find that xi — X2- So we conclude that 
the energy shift due to Stark effect is the same for any two magnetic field insensitive 
states. 



Appendix B. Driving stimulated Raman transitions using an Electro-optic 
Modulator 

The fields driving stimulated Raman transition in ions are typically generated from 
a single laser source with the multiple frequencies generated by optical modulators. 
Acousto-optic modulators can produce frequency shifts up to about 1 GHz, while 
electro-optic modulators can modulate at upwards of 10 GHz or more. The electro- 
optic modulator offers a solution for driving Raman transitions in ion species with a 
large hyperfine splitting, but unlike acousto-optic modulators, all frequencies in the 
modulated field have the same wave-vector, making it difficult to separate different 
frequency components. 
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Electro-optic modulators control the birefringence of uniaxial crystals with a lower 
frequency electric field, effectively modulating either the phase or the polarization of 
the incident optical field, depending on the orientation of the optical axes. Since the 
Raman coupling is polarization-dependent, the polarization modulation is equivalent 
to an amplitude modulation. For non-copropagating Raman beam geometry, the 
modulated field is divided using a beam splitter and the two beams recombine at the 
trap from different angles. If the difference between any frequency from one beam and 
any frequency from the other beam matches the energy splitting of two atomic and/or 
phonon levels, then a transition can potentially be driven. However, the amplitude 
of each pair of frequencies driving the transition can result in cancellations in the total 
transition rate. Usually, amplitude modulation produces sidebands with the same phase, 
resulting in the transition rates adding constructively. But phase modulation produces 
a comb of sidebands, some having amplitudes with opposite phases, which could result 
in a total transition rate of zero. This problem can be remedied by setting the beam 
path length difference between the two beams to certain values that produces a non- 
zero total transition rate [32] • This effectively adds a different phase to each sideband, 
resulting in a transition rate proportional to the squared electric field 

oo 

~ Yl U<P)e me Jn + k(<py in+k)d = J k (20 sin (9)) (B.l) 

n=— oo 

where J n (x) is the n-th order Bessel function. Equation IB. II describes the Raman 
transition rate involving the optical carrier and the fc-th frequency modulated sideband 
with modulation index <fi with a phase shift of 9 = (5kAx)mod(27r). Here 5k is the 
wave-number associated with the modulation frequency and Ax is the beam path length 
difference. 

To avoid nulling the average intensity due to destructive interference between the 
two fields at the ion, the field propagating along one path can be frequency shifted 
slightly from the field in the other path. The modulation frequency would then have to 
be adjusted to compensate for this frequency shift so that pairs of frequencies matches 
the energy difference between the two coupled levels. When the frequency offset is 
accounted for, all the Ak vectors driving the transition have the same sign, and the 
resulting transition is exactly the same as if each of the two beams having only a single 
frequency. To reverse the Ak vector, only the frequency offset on one beam needs 
to be changed. For example, to drive a Raman transition between two levels which 
have an energy splitting of h\j t ransition, the modulation frequency of the EOM can be 
set to ujeo < ^transition- The modulated beam is split into two paths, with the beam 
in path A frequency shifted by u ff set and the beam in path B frequency shifted by 
^offset + ^transition — ^eo- The two beams are recombined at the ion, with wave- vector 
\tA and kg respective to beam path A and B. In this case the wave- vector difference 
Ak for the Raman transition is equal to kg— k^, since the beam in path B has higher 
frequency. The beam in path B can also be frequency shifted by w ff set +uj E o — oj transition 
instead, in which case the wave- vector difference Ak would be k^ — kg since the beam 
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in path A would have the higher frequency. 

The reversal of Ak is useful in the phase stable configuration for a a$ gate 
(see section 2.3.2). To generate red and blue sidebands with opposite wave- vector 
difference Ak r = — Ak& , the frequency of the field along kg can be shifted by 
Uoffset+u' — co> 2 — 8— ojeo to generate the red sideband and by oj ff set —uj' —oj2—5+u)EO to 
generate the blue sideband. These two frequencies can be arbitrarily close to one another 
by tuning the modulation frequency of the EO close to the qubit frequency splitting uj' , 
which allows both frequencies to be generated using a single frequency shifter with a 
given bandwidth. However, if the modulation frequency of the EO is exactly u , then 
each beam would simultaneously drive a copropagating carrier transition (see section 
1.2.1) that will interfere with the operation. Therefore the modulation frequency 
should be tuned to approximately but not exactly u' . 
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Figure 1. 111 Cd + as an example of a hyperfinc qubit. The 2 S1/2 ground state 
electron configuration combined with the nuclear spin / = 1/2 creates an energy 
splitting of 14.5GHz due to the hyperfinc interaction. The states \F — 0, mp = 0) 
and \F = 1, mjr = 0} are called "clock" states at zero magnetic field since their energy 
difference has no first order dependence on the magnetic field, and are designated 
as the qubit states ||) and |!) respectively. A small external magnetic field lifts the 
degeneracy of the F = 1 states through the Zeeman effect. Due to the large ground 
state hyperfine splitting, the qubit can be initialized by optically pumping into the 
|t) = \F = 0,m,F = 0) state. The qubit state can be measured by applying resonant 
<j + radiation to optically pump the ||) = \F = 1, rap = 0) state to the \F — 1, mp = 1) 
state and drive a cycling transition between the 2 Si/ 2 , F = 1 and the 2 P^/2, F = 2 
excited state and collecting the resulting fluorescence. 
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Figure 2. Qubit with coherent Raman coupling via an excited state. Fields at 
frequency uu a and top couple the qubit levels |f) and \l) via the excited state |e). 
The fields are detuned from the excited state resonances u>^ e and Co^ e by frequency 
A. 
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Figure 3. Stimulated Raman transition between vibrational levels. Coupling depends 
on the beat note of the two Raman fields u>p — uj a : a) lo' for carrier transition, b) u'q—lu^ 



for first red sideband transition, and c) loL 



for first blue sideband transition. The 



qubit frequency splitting shifts from luq to u' due to AC Stark effect when the fields 
are turned on. 




Ion#l 



Ion #2 



Ion#l 



Figure 4. Cirac-Zoller Gate Scheme jj. A phase gate is constructed by performing 
the sequence illustrated here: i) a 7r-pulse on the first blue sideband on the first ion 
to map the internal state to the collective vibrational state; ii) a 27r-pulse between the 
|J,,n= 1) state and an auxiliary state \aux) on the second ion, resulting in a 7r phase 
shift on the state \l,n = 1); iii) a 7r-pulse on the first blue sideband on the first ion 
to map the vibrational state back to the internal state. A controlled-NOT gate can 
be constructed from a phase gate with a 7r/2-pulse on the second qubit before and 
after the phase gate. While the phase gate by definition maintains strict control of the 
qubit phase, the controlled-NOT gate relies on the two additional 7r/2-pulses having a 
particular phase with respect to the qubit. 




Figure 5. Representation of the eigenstate of a) a ■ z and b) <x • <fi on the Bloch sphere 
corresponding to the eigenvalue +1. A spin-dependent force creates two separate 
coherent states in phase space corresponding to the eigenstates of cr ■ n, as represented 
in c) , thus entangling the internal spin with the external motion of the ion. 




Figure 6. a) Displacement in phase space, in a frame rotating at the natural frequency 
of the harmonic oscillator. The displacement operator translates motional states in 
position/momentum phase space without distortion, b) For a force detuned from 
resonance, the motional state follows a circular path. For a closed trajectory, the 
quantum state acquires a geometric phase cj) = A/h in & round-trip orbit, where A is 
the area enclosed by the trajectory. 
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Figure 7. A er z -dependent force is driven by electromagnetic fields with two 
frequencies separated by u> v +8, as shown in a). These fields couple the two qubit states 
to the excited states with different coupling strengths (depending on polarization), 
producing a differential AC Stark shift that oscillates at u u + 6. The two fields must 
have a non-zero wave-vector difference Ak = k^ — k^ with a component in the x 
direction. 
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Figure 8. The spacing of the two ions determines the relative phase of the optical 
field experienced by each ion. a) The ions drawn are spaced by an integer multiple of 
the optical wavelength, creating an equal force on the two ions given the same internal 
state. In this scenario, b) spin states with opposite parity can excite the stretch mode, 
and the geometric phase acquired by each state is proportional to the area covered 
by the trajectory in phase space (tt/2 for the phase gate shown in Equation I27J) . The 
spin states with the same parity remain at the origin (not shown here) and acquire no 
geometric phase. 
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Figure 9. A er^-dependent force is driven by electromagnetic fields with at least three 
optical frequencies as shown in a). Two frequencies separated by lu' — lo v — 5 drive a 
detuned red sideband and a third frequency differs from one of them by lu' + lu u + 5 
to drive a detuned blue sideband, i) and ii) are two examples of possible frequency 
configurations. Some of the fields can have overlapping wave- vectors, but any pair of 
frequencies that drives a sideband must have a non-zero wave- vector difference with a 
component in the x direction. 
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Figure 10. Possible beam configurations for the gate, a) Phase sensitive 
configuration. The wave-vector difference for pairs of frequencies driving the red 
sideband and the blue sideband travel in the same direction Ak r = Akb, using the 
frequency configuration shown in Figure E| a) i) ) . A phase shift 8<j> in one beam 
path results in a phase shift in the spin of the entangled state, b) Phase insensitive 
configuration. The wave-vector difference for pairs of frequencies driving the red 
sideband and the blue sideband travel in the opposite direction Ak,. = Ak;,, using 
the frequency configuration shown in Figure a) ii) . A phase shift 8<f> in one beam 
path results in no net phase shift in the spin of the entangled state. 
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Figure 11. Demonstration of the phase sensitivity of the M0lmer-S0rensen (ct^) force 
applied to a single ion using different beam configurations, (a) Pulse sequence of a 
photon-echo experiment for testing optical phase sensitivity of the force. The echo 
pulses are carrier transitions driven by non-copropagating Raman beams propagating 
along the same wave-vectors as the force, and an applied phase shift 4>o is added 
using an acousto-optic modulator that controls the timing of the pulse. The o$ force is 
applied for sufficient time such that the two motional states corresponding to spin states 
IT^) and I J.^) have very little overlap at time t. A separate pulse in the other arm of the 
echo experiment cancels the residual AC Stark shift induced by the field driving the 
force (the ratio wo/A is significant enough to produce a non- negligible differential Stark 
shift between the two qubit states in this experiment) . (b) Probability of detecting 
I J,) vs. applied shift in the phase of the echo pulses <po using the phase sensitive 
configuration described in section 2.3.1 and Figure ITUh. The fringe contrast shows 
coherence between the phase <fi in the 0$ force and the phase </>o in the Raman carrier 
pulses (Probability should vary sinusoidally from to 0.5 when there is no decoherence) . 
(c) Same plot using the phase insensitive configuration described in section 2.3.2 and 
Figure If Ob . This time there is no coherence between the phase-insensitive o§ force 
and the phase-sensitive non-copropagating Raman carrier pulses. 
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Figure 12. Raman spectrum of two ions in the trap using a frequency comb generated 
by an electro-optic modulator (modulation frequency uj-eq— ujq — 1.5MHz in this case). 
The x-axis shows the frequency difference between the fields along the two beam paths. 
The carrier transition appears at ±1.5MHz (C), with corresponding first center-of-mass 
blue sideband transition at =f0.6MHz (B\) (uji/2-k = 2.1MHz) , first center-of-mass red 
sideband transition at ±3.6MHz (Ri), first stretch mode blue sideband at =p2.1MHz 
(B 2 ) (uj 2 /2tt = 3.6MHz), and first stretch mode red sideband at ±5.1MHz (R 2 ). 



